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RELATIONSHIP BETWEEN FOREST CLEARING AND BIOPHYSICAL 
FACTORS IN TROPICAL ENVIRONMENTS: IMPLICATIONS FOR THE 

DESIGN OF A FOREST CHANGE MONITORING APPROACH 

Steven A. Sader and Armond T. Joyce 
NASA/Earth Resources Laboratory 
National Space Technology Laboratories 
NSTL, Mississippi 39529 

ABSTRACT 

The scope of the work embodied in this paper was to analyze the relation- 
ship between forest clearing, biophysical factors (e.g., ecological zones, 
slope gradient, soils) and transportation network. The country of Costa Rica 
was selected as a study site because of the ecological diversity present with- 
in a relatively small area, and the availability of maps depicting biophysical 
information. 

The location of forested areas at four reference dates (1940, 1950, 1961, 
and 1977) as derived from aerial photography and Landsat MSS data was digi- 
tized and entered into a geographical ly-referenced data base. Ecological 

zones as portrayed by the Holdridge Life Zone Ecology System, and the location 
of roads and railways were also digitized from maps of the entire country as 
input to the data base. Information on slope gradient and soils was digitized 
from maps of a 21,000 square kilometer area, corresponding to the San Jose 
1:200,000 Quadrangle map. The total area of forest cleared over four decades 
as related to biophysical factors, was analyzed within the data base and de- 
forestation rates and trends were tabulated. The relationship between forest 
clearing and ecological zone and vhe influence of topography, soils, and 
transportation network are presented a. id discussed. Quantitative information 
on these relationships and data from other countries will form the basis for 
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predictive modeling of changes in the forest environment in tropical and sub- 
tropical regions. 


Introduction 

Tropical forests have great Importance with respect to the global energy 
balance, climate, geochemical cycle and various beneficial resources they pro- 
vide to mankind. Although considerable economic and human resources have been 
allocated to study the envi ronmental Impacts of forest clearing or "deforesta- 
tion" in tropical regions (U.S. Department of State, IS78 and 1380, World 
Bank, 1978, Zerbe et al., 1980 and Barney, '980), there are many uncertainties 
concerning actual rates and trends of tropical deforestation. Tropical de- 
forestation estimates differ by 100% or more (Browi and Lugo, 1930). 

A National Academy of Sciences report (Myers, 1978) has been cited in many 
U.S. Government reports as the most complete work on "tropical deforestation", 
although some authors dispute Myer's deforestation estimates and related con- 
clusions concerning environmental impacts (Sedjo and Clawson, 1983). A Food 
and Agricultural Organization (FA0) report (Lanly, 1982) sponsored by the UNEP 
Global Environmental Monitoring Program (GEMS), may have replaced the Myers 
report as the "best available data" concerning tropical deforestation. 
However, like the Myers report that preceded it, there Is no quantifiable 
basis to determine the reliability of the FA0 estimates. The FA0 report could 
serve as a baseline against which future tropical forest change can be mea- 
sured in the 76 countries targeted by the document, but a more standardized, 
quantifiable approach Is needed to derive current deforestation estimates. 

Several authors have demonstrated the utility of remotoly sensed data for 
detecting foresx change (Sader 1980; Joyce et al., 1980; Miller et al, 1979 
and Novaes, 1980). Remote sensing techniques coupled with computer technology 
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and appropriate sampling methodology may offer the best opportunity to system- 
atically monitor changes and trends in the tropic^ 1 (frost-free) environments. 

NASA's Remote Sensing Research in Tropical Environments 

A research program entitled “Remote Sensing of Forest Dynamics In Tropical 
Regions" is being conducted at the NASA-Earth Resources Laboratory of the 
National Space Technology Laboratories in Mississippi. The research addresses 
two general objectives: (1) to develop an ecologically-oriented stratifica- 

tion scheme and (2) to design a tropical forest inventory and monitoring 
system utilizing remotely sensed and biophysical (terrain, soils, etc.) data 
analyzed in a computerized Geographic Information System (GIS). The ultimate 
goal of the projected 4 year research program is to develop globally appli- 
cable sampling/inventory methcdologies to increase the precision of biorass 
estimates, timber volume, organic carbon and COo flux, and current deforesta- 
tion rates and trends. 


Approach 

The first task of NASA's investigation was to build a data base and 
analyze the relationship between forest clearing and biophysical landscape 
characteristics. The country of Costa Rica was selected as a study site 

because of the ecological diversity present within a relatively small area and 
the availability of published maps depicting biophysical information. The 
major hypothesis set forth in this investigation was that the propensity of a 
forest landscape to be cleared of primary forest is related to various bio- 
physical characteristics of the landscape and other man-induced factors. 
Specific hypotheses concerning forest change by life zone, soils and slope 
class will be investigated to address the major hypothesis stated above. 
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Furthermore, by quantifying these relationships and Identifying Important 
"predictors" of forest change, an Improved forest change Inventory and moni- 
toring approach can be Implemented utilizing remotely sensed (orbital satel- 
lites) and corrollary data. 

Methods 

Forest maps depicting the location and extent of all primary, dense, (>80% 
crown closure) forest In Costa Rica at four reference dates were acquired from 
the Costa Rican National Planning Agency (Perez and Prottl, 1978). The forest 
maps were derived through interpretation of aerial photos (1940. 1950. 1961) 
and image analysis of Landsat MSS data (1977). The forest boundaries of eacf 
date were transferred to a base map, digitized and reformatted to raster data 
In the ERL computer. 

Also ecological zones of the entire country as portrayed by the Holdridge 
Life Zone System (Holdridge, 1967) and transportation network (roads and rail- 
roads) were digitized from published map sources. Information on slope grad- 
ient, soils and geomorphology were acquired from the Instituto Geografico 
National. In addition, data were entered into the data base for a 21,000 
square kilometer area that was selected for detailed study, corresponding to 
the San Jose 1:200,000 quadrangle map sheet. 

All digitized data were transformed to raster image files in a common geo- 
referenced data base. In this form. Individual files could be queried to gen- 
erate tabular summaries of land area occupied by each class of Information. 
Manipulation and merging of the data files was performed to examine the rela- 
tionship of various biophysical and man-induced factors to the phenomena of 
forest clearing. 
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Results and Discussion 


At the country level, tabular summaries of the total extent of forest area 
were obtained at each reference date. These data Indicate that total primary 
forest remaining In 1940, 1950, 1961, and 1977 was 67X, 56t, 45X. and 32%, 
respectively. Assuming that Costa Rica was originally 99. 81 forested (Trop- 
ical Science Center, 1982), app, oxlmately 2/3 and 1/3 of the primary forest 
remained In 1940 and 1977, respectively. Average annual rates of deforesta- 
tion In primary forest Increased slightly from the 1940-50 period (1.6t/year) 
to the 1950-61 period (1. bi/year) and again In t'c 1961-77 period (1.9t/year). 
Specific relationships and trends In forest clearing will be examined in the 
following series of data analyses. 

Forest Clearing vs. Life Zone 

The Life Zone map of Costa Rica contained 20 life zones with 8 being tran- 
sitional from one zone to another. Each life zone is characterized by the 
Integration of 3 climatic parameters: potential evapotranspiration, precipi- 
tation and biotemperature (Holdridge, et al., 1971; Tosi, 1969). Forest 
clearing rates and trends by life zone (excluding Subalpine Rain) were 
examined for four time periods (pre 1940, 1940-50, 1950-61, and 1961-77). For 
all time periods except 1961-77, the largest amount of forest area was removed 
In the Tropical Moist Life Zone, followed by the Premontane Wet Zone (Table 
1). For the 1961-77 period, the Premontane Wet Forest was highest (2665 km 
removed) and the Tropical Wet was second highest In forest land cleared. 
Primary forest occupied only lit of the entire Tropical Dry Life Zone as of 
1940. Maps depicting forest area prior to 1940 were not available to allow 
estimates of forest clearing rates in the Tropical Dry and other Life Zones 
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for that period. The highest average annual rates 1 of forest clearing 
occurred In the T.oplcal Dry (9.9X) and the Premontane Moist Forest (7.91) 
between 1940 and 1950. The Premontane Moist and Tropical Moist Life Zones had 
the highest annual rates of deforestation In the last three time periods. 

Tosl (1980) reported that between 1971 and 1977, among Life Zones with a 
large total land area, the Tropical Mr 1 si Life Zone had the highest annual 
rate of deforestation (7.6X), followed by the Premontane Wet (3.7X) and the 
Tropical Wet (2X). Our computations (1961-77) show the same three life zones 
In the same order as having the highest deforestation rates, although our 
annual rates are slightly lower over the longer period. Tosi's estimates 
Include both primary and secondary forest where our computations relate only 
to clearing of primary forest. 

Forest clearing within the land area represented by the San Jose 1:200,000 
scale quadrangle show similar trends to the entire country results (Table 2). 
Of the Life Zones with a larger land area, the Tropical Moist Forest had the 
highest annual forest clearing rate for the 1940-50 and 1950-61 periods. The 
Tropical Moist was one of 4 life zones that were totally cleared of primary 
forest by 1961. The Tropical Wet and Premontane Wet life zones had the 
greatest amount of forest area cleared In the last 3 time periods. Annual 
forest clearing rates were slightly higher for the San Jose quadrangle as 
compared to the country level. For the 1950-6 1 period, the annual rate was 
2.7X for the San Jose quadrangle compared to 1.8X for the entire C "try. The 
San Jose quadrangle contains the majority of the urbanized area and popuTation 
within Cost* Rica. Population density by Life Zone according to Tosi and 

J The Montane Wet and the Tropical Dry Life Zone had a 9.1X/year rate of clear- 
ing but the total area cleared was only lOknr and 4km, respectively. 
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Voertman (1964) was highest in the Lower Montane Moist (238.5 persons/km*) . 
The Premontane Moist Life Zone had the next highest population density 
(216 .2/km 2 ) . Table 1 and 2 indicate total forest loss In the Lower Montane 
Moist and greater than 99J loss in the Premontane Moist Life Zone by 1961. 

Up to this point we have discussed the Life Zones that have experienced 
the highest forest clearing rates. There are some interesting trends in Life 
Zones that have experienced low amounts and rates of forest clearing that may 
be correlated with the life zone parameters. The hypothesis concerning these 
data are that forest clearing decreases as the potential evapotranspiration 
(PET) to precipitation (P) ratio decreases . Climates are increasingly humid 
as the ratio decreases. Although country-wide, the three life zones with the 
lowest PET/P ratio were the three that had experienced the lowest amount of 
forest clearing, the Montane Rain was the least deforested in 1977 (121 loss) 
but had a slightly higher PET/P ratio than both the lower Montane and 
Premontane Rain forest (Table 1). For the intensive study area corresponding 
to the San Jose 1:200,000 quadrangle map sheet, the hypothesis appears to be 
especially valid, although not tested in the rigid statistical sense in this 
preliminary examination. The Premontane Rain Forest was the least oefo^ested 
Life Zone in 1977, followed by the Lower Montane Rain and Montane Rain Forest. 
Likewise, the Premontane Rain Life Zone represents the lowest PET/P ratio of 
all life zones in Costa Rica followed by Lower Montane Rain and Montane Rain 
Forest, respectively (Table 2). 

A corollary hypothesis to be Investigated at a later stage In the research 
is that the probability of abandonment (and secondary succession) following 

2 

The reader is referred to Holdridge (1967), Holdridge, et. al. (1971), and 

Tosi (196?' for more information on life zone characteristics. 
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forest clearing will increase as the PET/P ratio decreases because the wetter, 
more humid zones are less suitable for non-forest uses and are prone to de- 
grade quickly. Brown and L igo (1980) have reported that mean plant biomass (in 
relatively undisturbed forest conditions) was linearly correlated to the PET/P 
ratio. Deforestation estimates by life zone would be a particularly useful 
statistic if biomass levels can be correlated with the PET/P ratio as Brown 
and Lugo suggest. 

Forest Clearing vs. Slope Gradient 

The digitized slope map contained six slope gradient classes. Computa- 
tions performed by computer indicated that 61% of the San Jose quadrangle was 
composed of greater than 30% slopes and only 5% of the slopes were nearly 
level. By 1940, the 30-45% slope class had the largest forest area cleared 
(Table 3). Between 1940 and 1950, the highest forest clearing rate occurred 
on 5-15% slopes (3.0%/year) and the lowest rate occurred on the 0-5% slope 
class (.4%/year). Between 1950 and 1961, the 0-5% slope class had only 
.6%/year clearing rate while the next lowest rate jumped to 2.3%/year for 
15-30% slopes. The rest of the slope classes were near 3%/year. Between 1961 
and 1977, the forest clearing rate dropped considerably on the steeper slope 
classes and increased substantially on the nearly level slopes (4.4%/year). 

The high rate of clearing on the steeper slope classes during the 1950-61 
period may be explained by the expansion of the cattle ranching industry 
during the same period (Tropical Science Center 1982, Tosi 1980). The de- 
crease in forest land cleared on the steeper slopes may have been a result of 
conservation awareness that materialized In the 1970's and the designation of 
new national parks, and forest protection zones (e.g., Braulio Carrillo, Rio 
Macho reserve and Zona Protectora South of La Selva Biological Reserve). 
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The total forest loss In relation to slope class by 1977 (Table 3), was gen- 
erally a linear decrease in the percentage of forest cleared as slope gradient 
increased. The trend was subtle up to the 60* and greater slope gradient 
class where an abrupt decrease occurred (e.g., 79* loss on 45-60* slopes 
compared to 54* loss on greater than 60* slopes). This relationship between 
forest clearing and slope gradient does not approach linearity until the last 
period. The steepest slope class (>60X) could be considered critical manage- 
ment areas; because the soils are prone to erosion, degradation and subsequent 
abandonment after the protective forest cover is removed. 

Forest Clearing vs. Soils 

Seven soil orders and 39 subgroups were present in the San Jose quad- 
rangle. The seven soil orders were: 1) Alfisols, 2) Entisols, 3) Histosols, 
4) Mollisols, 5) Ultisols, 6) Vertisols, and 7) Inceptisols. The Alfisols and 
Vertisols had been totally cleared of forest cover by 1940. By 1977, no 
primary forest remained on the Mollisols. The total area represented by 
Alfisols, Vertisols and Mollisols was small and they occur in the Tropical Dry 
and Moist Life Zones which were favored for agricultural and grazing uses. 

The largest amount of land area and some of the highest annuil ratas of 
forest clearing occurred on two soil subgroups (Typic Dystrandept-It ard Typic 
Dystropept-126) . They are also the two Icnjest in terms of total land area 
(Table 4). The lowest rates of clearing occurred on Typic Placandepts (14), 
Typic Hydraquents (El) and Oxic Palehumults (Ul) which had moderate to strong 
use limitations for cultivated crops and pasture as a result of either rocky, 
shallow soils, high rainfall, poor internal drainage or high soil acidity. 
The Typic Placandept was the only soil unit that experienced no loss in forest 
cover. It was learned by pinpointing the soil location that the reason for 
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Table 4. Forest Area by Soil Class at Four Itefprence Data - San Jose 1 : 200 .000 Quadrangle 
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(•Soil Association descriptions are by Peret, et al. 1479) 


it's unique status (zero deforestation) was probably not because there were 
inherent soil limitations for use, but because the mapping unit was totally 
enclosed by the boundaries of Rio Macho Forest Reserve. The significance of 
soils to forest clearing and the relationship of soils to other biophysical 
landscape characteristics will be Investigated in more detail as the research 
progresses into the next stage. 

Forest Clearing vs. Transportation Network 

The transportation map depicts the road and railroad system as of 1977 
when the map was revised. A hypothesis concerning these data are that forest 
clearing increases immediately adjacent to new road (or railroad) construction 
and out to an undefined distance from the nearest road or railroad access. 
The approximate dates and sequence of road construction in Costa Rica can be 
deduced by observing the 1977 transportation network with each date of forest 
location. Point A and B (Figure 1) show roads crossing forest in 1940 but the 
forest adjacent to the road is gone in 1950. It is likely that the section of 
road at point A and B was constructed after 1940 but before 1950. Using the 
same logic, the road at point C and D was constructed between 1961 and 1977. 
Without attempting to separate cause and effect, there appears to be a defi- 
nite relationship between transportation infrastructure and forest clearing as 
of 1977. 

This simple example demonstrates how historic Information (transportation 
Infrastructure development) can be derived from remotely sensed data if ade- 
quate records are not available in other forms. Our interest concerning 
these data is to determine the relationship between transportation infra- 
structure and the phenomena of forest clearing. Once additional maps can be 
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acquired depleting the transportation system at earlier dates, a distance-to- 
road spatial analysis criteria will be Investigated as a possible component of 
a "deforestation event probability" model. Also, as we have attempted to 
demonstrate, the terrain, climate and various demographic variables will have 
significant Influence on the amount of forest cleared adjacent to roads and 
railroads. 


Summary and Conclusions 

The relationship of biophysical landscape characteristics to the phenomena 
of forest clearing In Costa Rica was Investigated. Various published maps 
depicting forest area, life zone, soils, slope gradient, and transportation 
network were acquired and entered Into a common geographic information system 
(see Figure 2). Tabular summaries generated by combining the forest area 
files with biophysical and transportation data files revealed specific rates 
and trends of forest clearing over a 4 decade period. In the most simplistic 
description of these results, we learned that less forest land was cleared (as 
of 1977): 

1. In Life Zones with the lowest PET/P ratios, and in Life Zones with 
lower population densities. 

2. As the slope gradient Increased. 

3. At locations not adjacent to roads or railroads. 

4. On soils with strong inherent limitations for crop and pasture uses. 

5. In areas assigned forest protection status such as national parks and 
biological reserves. 

The results are not particulary surprising and most would appear to be 
Intuitively obvious. The value of the automated GIS as demonstrated here was 
the ease in which data files could be combined and compared to Isolate Import- 
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Six slope classes are repre- 
sented. Yellow hues are lower 
slope gradients and green hues 
are steeper. The slope image 
was combined with the 1977 
forest area image that still sup- 
ported primary forest in 1977. 
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Each life zone is identified by 
the integration of three climatic 
parameters. The wetter, more 
humid life zones (purple and 
dark green hues) had the lowest 
clearing rates and were still 
forested in 1977. 


In general, the least amount of 
deforestation by 1977 occurred 
on soils that had inherent limita- 
tions for non forest uses (e.g., 
crops and pastures). 


Soils 


Forested Soils-1977 



1940 forest is shown in blue 
and 1977 is shown in green. 
Roads and railroads are red. 
There was an obvious relation- 
ship between the transportation 
system and forest clearing. A 
distance-to-road spatial analysis 
technique is being investigated 
to quantify this relationship. 
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ant predictors of forest or land use change, quantify their relationships, and 
pinpoint their location on the ground. Many of the variables affecting forest 
change are Interrelated and others were not evaluated In this preliminary 
analysis. Because the accuracy of the forest area maps are unknown, we pre- 
sent the area of primary forest cleared and deforestation rates only as approx- 
imations for comparative purposes. Also, the deforestation rates had to be 
averaged over each reference period. Actual rates could have been higher or 
lower In Intermediate years depending on economic factors and population 
migrations. 

The Interactions of biophysical variables and their relationship to forest 
change will be further examined at the next stage In the research where "pre- 
dictors" of forest change will be determined quantitatively and applied In a 
regression model using statistical programs similar to those reported by 
Hamilton and Wendt (1975) and Hamilton (1974). A "deforestation event proba- 
bility" model will be based on empirical relationships between biophysical 
variables and other predictors of forest change as derived from the Costa Rica 
data and data bases generated for other tropical and subtropical areas. 

The deforestation model will be applied and tested in other tropical re- 
gions where the "sampling frame" can be generated through the acquisition of 
satellite-acquired data and good quality maps depicting the biophysical land- 
scape and other data necessary as input to the model. In an operational mode, 
forest change statistics compiled by forest type, life zone, slope gradient, 
and soils, would be accompanied by confidence limits. Critical areas (e.g., 
steep slopes, erodible soils, endangered habitat) and secondary forest succes- 
sion could be monitored systematically with long-term environmental changes 
and trends documented through the use of data collected by satellite-borne 
sensors to update the GIS. 
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